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ABSTRACT 



We examine the physical parameters that affect the accumulation of gas in molecular clouds to high column densities where the 
formation of stars takes place. In particular, we analyze the dense gas mass fraction (DGMF) in a set of self-gravitating, isothermal, 
magnetohydrodynamic turbulence simulations including sink particles to model star formation. We find that the simulations predict 
close to exponential DGMFs over the column density range A'CH?) = 3 - 25 x 10^' cm"^ that can be easily probed via, e.g., dust 
extinction measurements. The exponential slopes correlate with the type of turbulence driving and also with the star formation effi- 
ciency. They are almost uncorrelated with the sonic Mach number and magnetic-field strength. The slopes at early stages of cloud 
evolution are steeper than at the later stages. A comparison of these predictions with observations shows that only simulations with 
relatively non-compressive driving (b < 0.4) agree with the DGMFs of nearby molecular clouds. Massive infrared dark clouds can 
show DGMFs that are in agreement with more compressive driving. The DGMFs of molecular clouds can be significantly affected by 
how compressive the turbulence is on average. Variations in the level of compression can cause scatter to the DGMF slopes, and some 
variation is indeed necessary to explain the spread of the observed DGMF slopes. The observed DGMF slopes can also be affected by 
the clouds' star formation activities and statistical cloud-to-cloud variations. 
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1. Introduction 

Star formation is ultimately controlled by the processes that 
regulate the formation of density enhancements in molecular 
clouds. In our current picture, the density statistics of the in- 
terstellar medium are heavily affected by superso nic turbulence 
(for a review, see iHennebelle & Falgaronell2012b . The density 
statistics depend on characteristics such as the total turbulent and 
magii e tic energy (e.g., lPadoan et al. 1997a; Nordlund & Padqap 
11999"; Vazquez-Semadeni & Garcia 2001; K owal et alj 120071 
iMolina et al. 2012; Federrath & Klessen 20lll FK13 here- 
after), the driving mechanism of the turbu lence (e.g., 
iFederrath et"aLll2010bt iFederrath & Klessen"2 012l FK12. here- 
after) , the equation o f state (e.g., Passot & Va zquez-Semadeni 
1998; Gazol & Kim '201 3|), and the driving scale (e.g., 
iFischera & Dopita.2004; Brunt et al.li2009.) . Constraining these 
characteristics is fundamental for virtually all analytic star for- 
mation theories. 

We have previously employed near-infrared dust extinc- 
tion mapp ing in analyzing columri density statistics of molecu- 
lar clo uds dKainulainen et al.l2009ll201 lallbtlKainulainen & TanI 
1201 3L KT13 hereafter). This technique is sensitive and well- 
calibrated at low column densities, making it suitable to study 
the mass reservoirs of molecular clouds. Exploiting this advan- 
tage, we studied how the clouds gather gas to the regime where 
star formation occurs. We used an easily accessible characteris- 



tic to quantify this, namely the dense gas mass fractiorQ (DGMF, 
hereafter), defined as a function that gives the fraction of the 
cloud's mass above a column density value 



dM'(> A^) 



M(>N) 



(1) 



where M(> A^) is the mass above the column density A^ and Mtot 
is the total mass. The DGMF is linked to the probability den- 
sity function (PDF), p{N), of column densities, which gives the 
probability to have a column density between [A^, A^ + dN], via 



dM' 



p(N')dN'/ p{N')dN', 



(2) 



where [Mow, A^high] is the probed column density range. The rea- 
son for analyzing DGMFs instead of PDFs is simply the intuitive 
connection to the total mass reservoir of the cloud. Previously , 
DGMFs have been analyzed by, e.g., iKainulainen et alJ ( l2009h 
who showed that starless clouds contain much less dense gas 
than star-forming clouds and by Ladaetal. (2010) who used 
them to derive a star-formation threshold. 

From the theoretical point-of-yiew, the form of the DGMF 
can be controlled by any of the forces affecting the cloud's 
density structure. The key parameters describing these forces 
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' We purposefully use here the term "dense gas mass fraction" in- 
stead of "cumulative mass function" (CMF) from our previous works. 
This is to avoid confusion with the "core mass function" that is com- 
monly used in literature. 
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arqj i) the sonic Mach number, Ms, H) the turbulence driv- 
ing (iFederrath et al.ll2008l l2010b ?). which is commonly denoted 
by b, with b = 1/3 corresponding to purely solenoidal driv- 
ing and b = I to fully compressive driving, and Hi) the mag- 
netic field strength, B, reflected by the Alfven Mach num- 
ber. Ma- These parameters relate to density fluctuations via 
(iNordlund & Padoanll999llPrice et alJioTTtlMolina et al.l2012h 



yS 



), 



(3) 



where (Tinp/(p} is the standard deviation of logarithmic, mean- 
no rmahzed densities and /? = 2M\/M^. This form of Eq. 
[3] dMolina et alj 1201 2|) is valid up to moderate magnetic field 
strengths. Ma ^ 6. The strength of the Ms - density cou- 
pling is of great importance for analytic star formation theo- 
ries, because it directly affects the star formation rates and - 
efficiencies (SFE) they p redict (e.g.,[K rumholz & McKe j[2005t 
iHennebeUe & Chabried 120111 jPadoan & Nordlundi i201ll see 
FK12). 

In this work, we will estimate how the different physical pa- 
rameters affect the observed DGMFs of molecular clouds. To 
this goal, we will analyze numerical turbulence simulations and 
derive predictions for observable DGMFs. We will then compare 
the predictions to t he results of Kain ulainen et al., (|2009 . 2011b,) 
and KT13 (see also lLadaet ani2010h . 

2. Simulation data 

We analyze a set of magneto-hydrodynamic simulations of 
isothermal, driven turbulence in a periodic box, including self- 
gravity and sink particles to follow gas accretion onto proto- 
stars (see FK12). Each simulation is a time-series which starts 
(t - 0) when the turbulence is fully developed and the gravity 
is switched on. Then, the evolution is followed as a function of 
SFE, defined as the fraction of mass accreted into sink particles. 
The formation of the first sink particle occurs at SFE - 0%. The 
sink particles affect their surroundings because of gas accretion, 
and we eliminated them from the simulations. The issue is de- 
scribed in Appendix lA.ll Here we quote the main result: the 
DGMFs of Ms = 10 simulations (which we directly compare 
with observations) with 512^ cells are unaffected by sink parti- 
cles below A^(H2) < 11 X 10^' cm"^. They are 70% accurate up 
to A^(H2) ~ 25 X 10^^ cm"^. We also show in Appendix lA!2l that 
the resolution does not affect the DGMFs in this range. 

The simulations were scaled so that their virial parame- 
ters, ffvir.o = 5crlL/{6GM) where o-y is the 3-D velocity dis- 
persion and L the size of the simulation, were close to unity. 
Observations have sh own that molecul ar clouds, on average, 
show avir.o ~ 1 (e.g., i Hever et al]|20 09^. However, this defini- 
tion is an idealized approximation. The actual virial parameters, 
ttvir = 2|£']iin|/|£'pot|, Vary by more than an order of magnitude 
in the simulations. However, the actual virial parameters do not 
affect density PDFs greatly (FK12). If the virial parameter is 
"low-enough" to allow some collapse, the density structure is 
determined by other parameters (FK12; Molina et al. 2012). 

To make a realistic comparison with observations, we pro- 
cessed the simulations with Ms = 5 - 10 to mimic data derived 
using near-infrared dust extinction mapping (iLombardi & Alvesl 
1200 1). First, column density data from simulations was re- 
gridded to 60"/pixel and smoothed to the FWHM - 120" res- 
olution (0.09 pc at 150 pc distance). The native resolution of the 



^ However, see the discussion on the caveat related to the Reynolds 
numbers of simulations in SectionlBTTI 



simulations with A^s > 10 is coarser than this, and we could not 
smooth them (we do not compare them with the lower Ms sim- 
ulations). Then, the column densities outside A^(H2) = [3,25] x 
10^' cm"^ were discarded, approximating the dynamic range of 
extinction mapping. The lower limit of the range was chosen to 
be high enough that it is possible to define separate structures 
in simulations using (approximately) closed contours of con- 
stant column density. This is because observ ationally, "clouds " 
are commonly defined in this manner (e.g., iLadaet al]|2010l) . 
Finally, Gaussian noise with o-(N) = 0.018A^(H2) + 0.2 x lO^' 
cm~^ was added, follow ing typical uncertainties in the data of 
Kai nulainen et alj (12009 ). This procedure was repeated for three 
different projections of the simulation data, and the DGMFs 
from them were averaged to form the final DGMF 

We examined the effects of the resolution and noise to the 
DGMFs. We experimented with the resolution of 0.03 pc that 
studi es employing Hersche l data of nearby clouds will reach 
(e.g., [Schneider et al.l l2013h . Similar resolution is reached by 
combined near- and mid-infrared extinction mapping when ap- 
plied to infrared dark clouds (IRDCs, KT13). The effect of the 
resolution and noise to the DGMFs was practically negligible. 



3. Results and Discussion 

3.1. Dependence of the DGMF on physical parameters 

We derived the DGMFs for the simulations up to SFE = 10%. 
Figure [T] shows the DGMFs of four simulations with Ms = 10 
and b - {1/3,0.4, 1). For the case b - 0.4, a non-magnetized 
and magnetized simulation is shown. The DGMFs at early stages 
(f = and SFE = 0%) are well-described by exponential func- 
tions, dM' oc e"'^ . When star formation begins, the DGMFs flat- 
ten. Their shapes remain close to an exponential function, or 
curve upwards approaching a powerlaw shape. This behavior is 
similar in all models. Since the DGMFs are close to exponential 
functions in the range A^(H2) = 3 - 1 1 x 10^' cm"^, we quantified 
their shapes through fits of exponentials. This yielded the range 
a = [-0.41,-0.023] in afl models. 

We examined the dependence of the DGMF slopes on the 
driving of turbulence and magnetic field strength (E) in the sim- 
ulations with Ms - 10. The results are shown in Fig.|2](left and 
center). Most importantly, the DGMF slope responds most sensi- 
tively to the turbulence driving, changing by a factor of 4.8 - 8.5 
when b changes from 1 /3 to 1 . The slopes depend clearly less on 
B. The non-magnetic simulations show significantly shallower 
slopes than magnetized ones, but if B > 3 yuG, the slopes are 
uncorrelated with it. 

The DGMF slopes depend on the SFE. The dependency is 
stronger in magnetized than in non-magnetized simulations: the 
spreads of the slopes in the range SFE = [1, 10]% for these cases 
are 0.09 and 0.03, respectively. The mean difference in the slopes 
of non-magnetized and magnetized runs is 0.05. The early stages 
(f = 0, SFE = 0%) show clearly steeper slopes than the higher 
SFEs. We also examined the relationship between the DGMF 
slopes and Ats- For this, we derived the DGMFs in the native res- 
olution of the simulations (smoothing would greatly reduce the 
size of the low-Als runs). Therefore, the results should be com- 
pared to observations with caution. Figure |2] shows the DGMF 
slopes and Als in simulations with b - 1/3. The slopes are non- 
responsive to Als, except when Ms - 5. 

The DGMFs can vary also due to /) the random nature of tur- 
bulence ("cloud-to-cloud" variations) and //) projection effects. 
The former can be examined by comparing simulations that have 
the same input parameters, but different random number seeds 



b=l/3, MIO, B=a 
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Fig. 1. DGMFs of four simulations (black lines) with Als = 10, processed to mimic those observed with near-infrared dust extinction 
mapping technique. The solid lines show the DGMFs at f = and the dotted lines at time-steps SEE - {1,3, 10)%. T he panels also 
show with dashed lines the mean DGME of nearby starless clouds (blue) and of Taurus (.Kainulainen et alj|2009l red), and of a 
sample of IRDCs (KT13, green). 



(e.g., #12, 14, and 17, see Table lATTl i. Unfortunately, we only 
had three simulation pairs with varying random number seeds. 
The mean difference in the DGME slopes among these was 0.08 
at the early stages (f - 0, SEE - 0%). However, for time-steps 
SEE > 1 the mean difference was only 0.02. The projection ef- 
fects were studied by examining the standard deviation of the 
slopes derived for three different projections of all models. The 
mean standard deviation of the slopes in all models was 0.03. 

We note that the effective Reynolds numbers of our sim- 
ulations (<10'*) are lower than that of the interstellar medium 
(~10^). It is not cle a r how this affects the predicted statisti- 
cal properties. lAluiel (1201 3l) has rigorously shown that the di- 
rect influence of driving on the kinetic energy is restricted to 
scales larger than the smallest scale at which the turbulence is 
stirred. However, numer ical (Eederrath et al. 2010b) and analytic 
dGaltier & Banerieg|201 II) works have found differences in flow 
statistics in the range that can be considered to be the "iner- 
tial range" of compressible turbulence simulations. Resolution 
studies of the simulations suggest that the driving-induced dif- 
ferences remain when the Reynolds number increases. As this 
issue cannot be addressed with the current computational meth- 
ods, our results are also subject to it. 



3.2. Comparing the predictions witii observations 

Figure s[T] and |2] show observed DGMEs to be compared with the 
simulated ones. Figure [1] shows the mean DGME of quiescent 
clouds (LDN1719, Lupus V, Cha III, and Mu sca) and a DGME 
of a typical star-forming cloud (Taurus) from iKainulainen et al.l 
(l2009i) . and a mean DGME of ten IRDCs from KT13. Figure [3 
shows the ranges of the observed slopes from Kainulainen et akj 
(120091) . which span a = [-0.17,-0.45] for 13 nearby star- 
forming clouds and a - [-0.35, -1.2] for four quiescent clouds. 
The range of IRDC slopes from KT13 is also shown. We note 
that the DGMEs of IRDCs in KT13 were derived from a slightly 
different column density range than those of nearby clouds (they 
begin from A^(H2) ~ 1 x 10^' cm"-). Thus, the comparison of 
them with the other data should be considered only suggestive. 

The dependence of the DGME slopes on the turbulence driv- 
ing allows us to constrain b (see Eig. |2]i. None of the simula- 
tions shows as steep slopes as observed in starless clouds. Erom 
the non-magnetized simulations, only those with b = 1/3 are in 
agreement with the nearby star-forming clouds. Magnetic fields 
can steepen the slopes by about 0.05 (Eig.|2] center). Therefore, 
from the magnetized runs those with b = 1/3, or fe = 0.4 and 
B > 3 fiG agree with star-forming clouds. The fully compres- 



sive simulations produce a greatly higher fraction of dense gas 
than observed in nearby clouds. The comparison suggests a low 
b for nearby molecular clouds on average, possibly lower than 
previously estimated bv Padoan et al. (1997b) and Brunt ( 20 lOl) 
in Taurus, b ~ 0.5. 

The DGME slopes correlate with the SEE, depending on 
whether the cloud is magnetized or not. Since in the current 
view clouds have magnetic fields (Crutcher 2012), the spread 
of slopes is likely the most realistic in magnetized simulations 
(i.e., 0.1, see Eig.|2]l. Thus, it seems that part of the spread in the 
observed slopes originates from the SEEs of the clouds. We used 
a Monte Carlo simulation to estimate whether all the variation in 
the observed slopes can originate from changes in the SEE and 
statistical variations. We assumed that the changes due to SEE 
are uniformly distributed between [0,0.1] and the statistical vari- 
ations are normally distributed with cr = 0.04. The test showed 
that the probability that 13 clouds span a range > 0.28 is 0.2%. 
Note that the range of the observed slopes can be wider KT13 
showed that IRDCs possibly have flatter DGMEs than nearby 
clouds (Eig.|2]i. In conclusion, it seems likely that the spread of 
the observed DGME slopes cannot be explained by statistical 
variations and changes in the SEE alone. Changes in the clouds' 
average compression provides one possible source to account for 
this variation. 

One interesting question for the future is to examine the ef- 
fect of cloud mass to the DGMEs. There are no very massive 
clouds in the nearby cloud sample (median mass 0.5 x 10'* M©). 
In contrast, the median mass of the IRDCs is 5 x 10"* M©, which 
is ten times higher This could contribute to the differences seen 
in the slopes of the two cloud sets. However as discussed earlier, 
comparing DGMEs of IRDCs with nearby clouds is not without 
caveats. The question could be properly addressed by a study of 
a statistical sample of IRDCs, or a study of the nearest high-mass 
clouds (e.g., Orion, Cygnus, Rosette) employing Herschel data. 

The weak dependence of the DGME slopes on Als appears 
to be an effect of the narrow column density range we examine 
(note that the results were derived from simulations that have dif- 
fering physical resolutions and are only suggestive). The density 
PDE is expected to respond to Als following Eq.|3] which should 
reflect to the DGMEs. However, it appears that in the range of 
A^(H2) - 3 - 11 X 10^' c m~^ the effect i s insignificant. This 
result is in agreement with 'Goo dman et aP (1200 9) who did not 
detect any dependence between column density PDE widths and 
CO linewidths in Perseus. However, we recently measured the 
column density PDE widths using a high-dynamic -range tech- 
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Fig. 2. Exponential slopes of the DGMFs as a function of h {left), B {center), and Als {right). The solid black lines show the timestep 
t - Q, and the dotted lines SFE - {0, 1, ... , 10}%. The b lue, red, and green shaded regions indicate the slopes observed in starless 
and star-forming nearby clouds (JKainulainen et alJ2009l) and in IRDCs (KT13, however, see the discussion on these data in Section 
I3.2l i. The median masses of each set of the clouds, M1/2, are shown in the panels. 



nique (KT13) and concluded that if a wider range is examined, 
the PDF widths correlate with Als- 

When comparing observed DGMFs with simulations, it 
should be kept in mind that in simulations "driving" is well- 
defined and ideal: energy is injected at large scales, with certain 
characteristics such as the divergence and curl. In real clouds, en- 
ergy is likely injected at multiple scales and the characteristics of 
the driving can depend on the scale. However, if some of these 
driving modes excite more compression than others, particular 
regions in a cloud, and hence, also clouds on average, can show 
characteristics of the flows produced with ideal driving with dif- 
ferent mixtures of solenoidal and compressive modes. 

Finally, we comment on the relation between the DGMFs 
and column density PDFs. The column density PDFs of nearby 
clouds are log-noimal below A^(H2) < 3x10^' cm"^. In the range 
A^(H2) = 3 - 25 X 10^' cm"^, they are in agreement with either 
powerlaws or (wide) log-normals (KT13). It is not established if 
the PDFs abov e MH?.) > 3 x 10^' cin '^ are log-normals (KT13) 
or powerlaws (ISchneider et al.ll2013[ see Fig. IB.U . Importantly, 
it follows from Eq. |2] that a log-normal PDF yields an expo- 
nential DGMF and a powerlaw PDF yields a powerlaw DGMF. 
The simulated DGMFs in the range A^(H2) > 3 - 25 x 10^' 
cm"^ appear exponential at the early stages. Therefore, the col- 
umn density PDFs at these stages are close to log-normals. When 
the simulations evolve, the DGMFs become closer to powerlaws 
(see FK13). This means that the underlying column density PDF 
transits from a log-normal to a powerlaw. 



4. Conclusions 

We have examined the relationship between the dense gas mass 
fraction (DGMF), star formation, and turbulence properties in 
molecular clouds by comparing DGMFs derived from isother- 
mal, magneto-hydrodynamic, self-gravitating turbulence simu- 
lations to observed ones. Our conclusions are as follows. 

1 . Simulations predict close-to exponential DGMFs for molec- 
ular clouds in the column density range of A^(H2) = 3 - 
11 X 10^' cm"^. The DGMF slopes span the range a - 
[-0.41, -0.023], being clearly steeper at the early stages of 
the simulations compared to the stages when stars are form- 
ing (SFE > 1%). These predictions are accurate on a 70% 
level up to A^(H2) ~ 25 x lO^' cm^^. 

2. The DGMF slopes depend strongly on the turbulence driving 
{b). They depend less, but significantly, on the exact SFE. 
The dependence on the SFE is stronger in magnetized than 



non-magnetized cases. Generally, the eflfect of the magnetic 
field to the DGMF is small. Also Als has a negligible effect 
on the slopes in the examined column density range. The sta- 
tistical variations are comparable to those arising from vary- 
ing SFE. However, how compressive the turbulence is (i.e., 
parameter b) is the largest single factor in determining the 
slope of the DGMF in the simulations. 
3. The obsei-ved DGMFs can be used to constrain the turbu- 
lence driving parameter b. The DGMFs of nearby clouds 
are only reproduced by simulations that are driven by rela- 
tively non-compressive forcing, i.e., b - 1 /3 or 0.4. The fully 
compressive simulations {b - \) over-estimate the DGMFs 
greatly. Massive IRDCs can show flatter DGMFs that are in 
agreement with more compressive driving. The spread of the 
obsei-ved DGMFs cannot be explained by different SFEs and 
statistical variations alone. Variations in the clouds' average 
compression level offer one explanation to account for the 
observed spread. 
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A. 1 . Effect of sink particles 

Sink particles (iFederrath et al.ll2010al) in the simulations accrete 
material into them after their creation, and hence, affect the den- 
sity structure of their immediate surroundings in the simulation 
(and the DGMFs). In the following, we consider the effects of 
sink particles to the DGMFs. 

As described in FK12, the sink particles are created on a 
certain, resolution-dependent volume density and always have a 
radius of 2.5 pixels in the native resolution of the simulation. It 
follows that the sink particles have a resolution-dependent min- 
imum density, which can further be converted into a minimum 
mean column density. Sink particles are created when a series 
of collapse criteria are fulfilled (see FK12), and when the local 
volume density exceeds 



Psink 



2 
Sink 



(A.l) 



where c^ is the isothermal speed of sound and rsink the radius of 
the sink particle. It follows that the mean column density of a 
sink particle at the moment of its creation is 



y _ Psink^sink _ 4 

■^sink — 9 — ^ Psink ^sink- 



(A.2) 



sink 



The sink particle properties are listed in Table lA. 1 I for different 
physical resolutions. 

The sink particle column densities listed in Table |AT| rep- 
resent levels below which the DGMFs are not affected by sink 
particles, regardless of whether the sinks are removed or not. 
In the most conservative interpretation, the DGMFs are reliable 
only below these column density limits. Therefore, we use the 
upper limit of A^(H2) = 11 x 10^' cm"^, which is the sink par- 
ticle column density for the A^s = 10 simulations 512^ cells in 
size, in the analysis performed in this paper 

However, it is not at all certain that the DGMF shape imme- 
diately above A^(H2)sink is greatly affected by the sink particles. 
Above A^(H2)sink, there are lines-of-sight whose column density 
is higher than the sink particle column density, but the local vol- 
ume densities do not reach high enough values for sink particles 
to form. In fact, these lines-of-sight are greatly more numerous 
in the simulations compared to those that contain sinks, espe- 
cially at early times when the overall SFE is low. 

We dealt with sink particles in this work by disregarding 
the lines-of-sight affected by them directly from the simulation 
data. While this procedure, in principle, eliminates the effects of 
sink particles, it removes mass preferentially from high column 
densities, and hence potentially biases the DGMF downwards 
(steepens it). Consequently, it is important to note that t\\t flat- 
tening of the DGMFs seen in the simulations (see Section 13.11 ) 
at around A^(H2) ~ 10-15x10^' cm"^ cannot be due to sink 
particle treatment; any associated incompleteness would bias the 
determination downwards, not upwards. 

We can quantify the incompleteness due to sink particle re- 
moval by comparing DGMFs derived with and without the elim- 
ination of sink particles. This experiment is shown in Fig. lA.ll 
which shows the ratio of the DGMFs with and without the sink 
particle elimination as a function of column density. The plot is 
shown for the model in which the effect of sinks in the exam- 
ined column density range is expected to be strongest, i.e., the 
solenoidal simulation with 256^ cell resolution. Higher resolu- 
tion increases the sink particle column density (cf.. Table IATTT i. 
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Table A.l. Simulation properties (adapted from FK13) 



# 


name|^ 


Ms 


b 


box size 


box size 


M 

m 


, 


ficL. 


Q'vir.O 


Q-vir 


n(H2).si„k 


A'(H2),i„k 


Msink 










[cells] 


[pc] 


[Mo] 


b"G] 






[10-* cm-3] 


[1021 cm-2] 


[Mo] 


6 


GT256sM5 


5 


1/3 


256 


2 


3.9 X 


W 





1.0 


8.0 


28 


22 


0.60 


7 


GT256mM5 


5 


0.4 


256 


2 


3.9 X 


10^ 





0.98 


5.4 


28 


22 


0.60 


8 


GT256cM5 


5 


1 


256 


2 


3.9 X 


102 





0.82 


1.5 


28 


22 


0.60 


10 


GT256sM10 


10 


1/3 


256 


8 


6.2 X 


10^ 





1.1 


12. 


1.7 


5.6 


2.4 


11 


GT512sM10 


10 


1/3 


512 


8 


6.2 X 


103 





1.1 


12. 


6.9 


11 


1.2 


12 


GT512mM10(sl) 


10 


0.4 


512 


8 


6.2 X 


103 





1.1 


4.5 


6.9 


11 


1.2 


13 


GT512mM10Bl(sl) 


10 


0.4 


512 


8 


6.2 X 


10' 


1 


1.1 


5.4 


6.9 


11 


1.2 


14 


GT512mM10(s2) 


10 


0.4 


512 


8 


6.2 X 


103 





1.2 


8.4 


6.9 


11 


1.2 


15 


GT512mM10Bl(s2) 


10 


0.4 


512 


8 


6.2 X 


103 


1 


1.2 


9.5 


6.9 


11 


1.2 


16 


GT256mM10(s3) 


10 


0.4 


256 


8 


6.2 X 


103 





1.0 


5.9 


1.7 


5.6 


2.4 


17 


GT512mM10(s3) 


10 


0.4 


512 


8 


6.2 X 


103 





1.0 


5.9 


6.9 


11 


1.2 


18 


GT512mM10Bl(s3) 


10 


0.4 


512 


8 


6.2 X 


103 


1 


0.97 


6.4 


6.9 


11 


1.2 


19 


GT256mM10B3(s3) 


10 


0.4 


256 


8 


6.2 X 


103 


3 


0.81 


8.4 


1.7 


5.6 


2.4 


20 


GT512mM10B3(s3) 


10 


0.4 


512 


8 


6.2 X 


103 


3 


0.83 


8.7 


6.9 


11 


1.2 


21 


GT256mM10B10(s3) 


10 


0.4 


256 


8 


6.2 X 


103 


10 


0.79 


6.6 


1.7 


5.6 


2.4 


23 


GT256cM10 


10 


1 


256 


8 


6.2 X 


103 





0.85 


1.1 


1.7 


5.6 


2.4 


24 


GT512cM10 


10 


1 


512 


8 


6.2 X 


103 





0.87 


1.1 


6.9 


11 


1.2 


25 


GT256sM20 


20 


1/3 


256 


32 


9.9 X 


10^ 





1.0 


11.1 


0.11 


1.4 


9.6 


26 


GT256mM20 


20 


0.4 


256 


32 


9.9 X 


10* 





1.1 


4.5 


0.11 


1.4 


9.6 


27 


GT256cM20 


20 


1 


256 


32 


9.9 X 


10* 





1.0 


0.60 


0.11 


1.4 


9.6 


28 


GT256sM50 


50 


1/3 


256 


200 


3.9 X 


10^ 





1.1 


12 


2.8 X 10-' 


0.22 


60 


29 


GT512sM50 


50 


1/3 


512 


200 


3.9 X 


10^ 





1.1 


13 


1.1 X 10-2 


0.44 


30 


30 


GT256mM50 


50 


0.4 


256 


200 


3.9 X 


10^ 





1.0 


7.0 


2.8 X 10-3 


0.22 


60 


31 


GT512mM50 


50 


0.4 


512 


200 


3.9 X 


10^ 





1.1 


7.4 


1.1 X 10-2 


0.44 


30 


32 


GT256cM50 


50 


1 


256 


200 


3.9 X 


10^ 





0.95 


0.54 


2.8 X 10-3 


0.22 


60 


33 


GT512cM50 


50 


1 


512 


200 


3.9 X 


10^ 





0.99 


0.56 


1.1 X 10-2 


0.44 


30 



" Parentheses after the names refer to the different random seeds used in the simulations. 

* Total mass in the simulation box. 

'' Mean magnetic-field strength in z-direction of the computational domain. 



and more compressive forcing increases the relative amount of 
high column densities, thereby reducing the error in the exam- 
ined column density regime. The figure shows that the error 
due to incompleteness (i.e., preferential removal of high-column 
densities) is less than 30% below A^(H2) < 25 x 10^' cm-^ for 
SFEs up to 10%. 

In summary, it can be concluded that the DGMFs derived for 
Ats = 10 simulations are unaffected by the sink particles (or by 
their removal) below the A^(H2).sink values. In addition, the error 
in the predicted DGMFs is less than 30% when the range up to 
A^(H2) Si 25 X 10^' cm-^ is considered. 

A.2. Effect of the simulation resolution 

The simulations of FK12 are either 128^ 256^ 512^ or 1024^ 
computational cells in size. In this work, we used all but those 
simulations that are 128^ cells in size. It is possible that the dif- 
ferent computational resolutions used in the simulations affect 
the DGMFs, as especially high column densities are potentially 
better resolved by higer-resolution simulations. We examined 
the possible effect of the simulation resolution to the DGMFs 
by comparing the DGMFs of simulations that were run with the 
same physical parameters, but have different computational res- 
olution. 



Figure IA.2I shows as an example a comparison of DGMFs 
derived for models #10 and #11 that are 256^ and 512^ cells 
in size, respectively. All other parameters are same in these 
two models. The figure shows the DGMF of the model #10 di- 
vided by that of model #1 1 (red line). The figure also shows 




° 0.70 



N(H ) [x 10^' cm- 



Fig. A.l. Error (incompleteness) in the derived DGMFs due to 
removal of sink particles. The figure shows the ratio of DGMFs 
derived with and without sink particle removal as a function of 
column density for time steps up to SFE = 10%. The curves 
for t - Q and SFE = 0% are indistinguishable from unity. The 
plot is shown for simulation #10 (Als = 10, 256^ cells in size, 
b = 1/3). The error in other Als - 10 models is expected to be 
smaller, because of the higher sink particle column density and 
more compressive turbulence driving. 



the DGMFs calculated for model #1 1 using different projections 
(projections to xy, xz, and yz planes, black dotted lines). The 
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Column density PDFs 



1.4 



^11 (512), (xiy, xz, yz) projections 
#10 (256), m|ean 




Fig. A.2. Effect of simulation resolution to the DGMFs. The red 
line shows the DGMF of simulation #10 (256^ cells in size) di- 
vided by the DGMF of simulation #1 1 (512^ cells in size). The 
physical parameters of the two simulations are the same. The 
dashed lines show the DGMFs calculated for different projec- 
tions of model #11 divided by the mean DGMF of model #11. 



DGMF of model #10 is in good agreement with that of model 
#11 below the sink particle column density, A^(H2) = 11 x 10^' 
cm"^. At higher column densities, the lower-resolution simula- 
tion (#10) begins to under-estimate the column densities slightly. 
However, it is still within 30% of the higher-resolution one up to 
the column density of A^(H2) ~ 25 x 10^' cm"^. We conclude that 
the effect of resolution is smaller than the uncertainty due to the 
projection effects in the column density range A^(H2) - 11x10^' 
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Fig.B.l. Column density PDFs of models #11 {b = 1/3) and 
#24 (b = 1), and the PDF of the Taurus molecular cloud. Both 
models have Ms = 10 and B = /iG, and they are 512^ com- 
putational cells in size. The black histograms show the PDFs of 
model #11 at f = (solid line) and SFE = 5% (dotted line). 
The blue line shows the PDF of model #24. The red line shows 
the PDF of Taurus from lKainulainen et aP (12009.) . Note that the 
dynamic range of the Taurus PDF ends at about In A^(H2) = 3.2. 
The black dashed line shows, for reference, a log-noimal func- 
tion. The PDFs in the range A^(H2) = 3 - 1 1 x 10^' cm"^ can be 
described by a log-normal function, but also reasonably well by 
a powerlaw function (which would be a linear curve in the given 
presentation). 



cm and accurate to 70% level up to A^(H2) 



25 x 10^1 cm-2. 



Appendix B: Illustration of column density PDFs 

Figure IBTI show a comparison of the column density PDFs de- 
rived for models #1 1 and #24, and t he PD F of the Taurus molec- 
ular cloud from iKainulainen et al.l (l2009l) . Note how the higher 
relative amount of high-column density material predicted by 
fully compressive simulations (#24) is evidenced by a flatter 
PDF In the column density range A^(H2) = 3 - 25 x 10^' cm"^, 
the PDF of simulation #1 1 is close to what is observed in Taurus. 
In this narrow range, the PDF is in a reasonable agreement with 
either a log-normal function (shown for a reference in the figure) 
or a powerlaw function. 



